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Abstract, in this talk I discuss three recent developments in the theoretical understanding of the phase diagram 
of the strong interaction. The first topic deals with the comparison of model calculations of the quark-hadron 
transition at vanishing quark chemical potential with state-of-the-art lattice QCD results. In the second relates to 
the size of a possible 'quarkyonic phase'. The third deals with the occurence of inhomogeneous chiral phases. 



1 Introduction 

Exploring the chiral and deconfining properties of strong- 
interaction matter at high temperatures and large densities 
is one of the central themes in nuclear- and astrophysics 
fl 2|. The general phase structure and, in particular, the 
possible existence of a (chiral) critical endpoint (CEP) for 
finite net baryon number density and its consequences for 
the phase structure of QCD at lower temperatures are much 
under debate. High-energy heavy-ion collision experiments 
at RHIC and the SPS have started to look for experimental 
evidence of the CEP and experiments at future heavy-ion 
facilities (FAIR and NIC A), have been designed to probe 
the relevant high-density region in the QCD phase dia- 
gram. 

Gaining insight into the properties of QCD matter at 
non- vanishing quark chemical potential (ju) is thus of great 
interest. However, the relevant region of the phase diagram 
is not easy to access directly in QCD. At non-zero values 
of jd the notorious fermion sign problem prohibits straight- 
forward lattice QCD simulations. Thus one currently re- 
sorts to models that capture some of the essential features 
of QCD, such as the basic symmetries and their break- 
ing patterns. Although limited in their predictive power, 
such models, nontheless, give valuable insight into various 
physical efiTects governing the phase structure, in particu- 
lar the existence of a chiral CEP, a confined but chirally 
restored phase at low temperatures ('quarkyonic matter') 
and/or the occurence of inhomogeneous chiral phases. 

Model predictions for the equation of state (EoS) of 
QCD matter can be compared with lattice simulations at 
vanishing ji. Due to large lattice sizes, small quark masses 
and algorithmic improvements, the latter have become very 
accurate recently. Such comparisons will be discussed in 
Sect. 2. Based on large Nc arguments for the T- and yu- 
dependence of the pressure, a 'quarkyonic' phase has been 
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conjectured at low T and large jn. In this phase, quarks and 
gluons are still confined in hadrons (mostly baryons) but 
chiral symmetry is restored, leading to parity degeneracy 
in the excitation spectrum. There has been much debate 
recently about the size and the detailed structure of this 
phase. Theoretical arguments and empirical constraints will 
be discussed in Sect. 3. An exciting possibility is the oc- 
curence of inhomogeneous chiral phases at low T and large 
yu, possibly entangled with inhomogeneous color supercon- 
ducting phases. Inhomogeneous phases are predicted in 
1 -h 1 dimensional models of QCD in the large N limit and 
their three dimensional analogue could shed light on the 
nature of the quarkyonic phase. Recent results and prelim- 
inary conclusions will be discussed in Sect. 4. 



2 The Quark-Hadron Transition at small fi 

I begin with the behavior of strong-interaction matter at 
very small yu and large T, as it is encountered in the early 
universe. Modern models such as the PNJL or the PQM 
model incorporate (approximate) chiral symmerty and its 
spontaneous breaking in the vacuum, the axial Ua(l) ano- 
maly and the heavy quark limit in terms of the temporal 
Polyakov loop. 



2.1 The Polyakov Quark Meson Model 

As an example I will discuss the EoS predictions of the 
PQM model (the results of the PNJL model are very simi- 
lar). The 3-flavor PQM model is basically the linear sigma 
model with up, down and strange quarks coupled the flavor 
octet of scalar and pseudoscalar meson fields cr^ and tt^: 

■Cquark = q iw - ((Ta + /yS^a)) ^ (1) 
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and the mesonic Lagrangian given by 

-Lmeson = Tv^d^M^d^M) - n?Tv{M^ M) + Tr[if (M + M^)] 
-^i[Tr(M^M)]2 - ^2Tr(M^M)2 + c (det(M) + det(M^)) 

(2) 



where 



M = ^^(cr, + /7rJ; if = ^ ^/z, . (3) 



The term involving the determinant of M and incor- 
porates the axial anomaly. This is supplemented by the 
Polyakov loop expectation value 



I = — TrPexp 



Jo 



dTA4(x,T) 



(4) 



and the covariant derivative of the temporal gauge field 
A4(x, r) as well as a potential term such that 

£poi = -q74A4q-m^J) • (5) 
The total PQM Lagrangian then reads 

-CpQM — -Cquark + -Cmeson + -Lpol • (6) 

Using the mean-field approximation, the parameters are 
adjusted to the vacuum meson masses and their weak de- 
cay constants. The Polyakov potential 14 {i, i), on the other 
hand, is determined from a fit to the pure gauge lattice EoS. 
There are various choices for 14 in the literature such as a 
polynomial in I and i based on Landau- Ginzburg theory 
O, a logarithmic ansatz motivated by the Haar measure 
of the SU(?>) color gauge group |4| or a choice derived 
from the strong-coupling expansion of QCD, proposed by 
K. Fukushima 13. 



2.2 Mean-Field Results 

At finite T and yu one evaluates the mean-field grand poten- 
tial as 

(7) 

where U (cri, cr^) denotes the eff'ective meson potential in 
terms of the light quark condensate cr/ = (//) and the strange 
quark condensate cr^ = (ss). The fermionic part Qqq in- 
volves both quark condensates and the Polyakov loop ex- 
pectation value and is given by 



Qqq = -INfT 



//) + ln^(r,-//)} (8) 



with 



where Ep = -h denotes the quark quasi-particle 

energy with constituent mass Mq. Eq. [9] shows that in the 
hadronic phase single- diquark contibutions are suppressed, 
since in the confined region both £ and £ vanish. Finally the 
phase diagram is determined via the stationarity condition: 



dQ 
dcTi 



dQ _dQ _ dQ 



= 0. 



(10) 



In Fig. [T] predictions of a PQM mean field calculation O 
are compared with recent lattice data of the HOTQCD and 
the WB colaborations 17181 . While the general agreement 
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Fig. 1. The normalized quark condensate A^s (upper part) and 
the Polyakov loop expectation value (lower part) as a function of 
T IT^ at vanishing ji. Here denotes the (pseudo)critical chiral 
transition temperature. The lattice data are taken from Ref. II7I8I 



is quite good, the quark-hadron transition is predicted to 
be sharper than in the lattice simulations (especially for 
physical bare quark masses |8 1). The same also holds for 
the pressure and other bulk thermodynamic quantities. We 
will return to this point in the next section. 



2.3 Including Fluctuations 

The relative sharpness of the quark-hadron transition in 
mean-field theory can be understood as an eff'ect of the 
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omission of quantum fluctuations. These can be included 
most efficiently in the functional renormalization group ap- 
proach (FRG). The FRG is based on an infrared regulariza- 
tion of the grand potential, which then becomes a momen- 
tum scale dependent quantity Qk with a scale parameter 
k. The full potential, including all quantum fluctuations is 
obtained in the limit /: ^ 0, i.e. 



Q(T,ii) = \imQk(T,/i) . 



(11) 



The momentum flow equations for the QM model have 
been derived in | 9 | and extended to include the Polaykov 
loop effects as a background field in [lOJ . The resulting 
ffow equations for Nf = 2 read 



dkOk 



127r2 

NcNf 



{\-n,{iJ)-n-,{iJ)) 



(12) 



Here Ub are the bosonic distribution functions with the pion 
and sigma energy 

eI = \+ 2Q[lk^ El = l-\- 2Q[lk^ + A(P^Q';ik^ (13) 

where the prime denotes a functional derivative with re- 
spect to the chiral condensate = (cr), i.e. 



Q[ = dQk/d(f) etc (p = {cr} . 



(14) 



The explicit form of the fermionic distribution functions in 
the Polyakov loop background field nq(£, i) and riqii, £) are 
given in 1 101 and involve the quark quasi-particle energy 



1 -h 



(15) 



In the absence of the background gluon field they reduce 
to the Fermi-Dirac distribution for constituent quarks, i.e. 
expression ([12]) coincides with that in (91. As can be seen 
from Fig.[2]the pure mean-field transition (left part) is sig- 
nifically softened when ffuctuations are included via the 
FRG. 
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Fig. 2. A comparison between mean-field results of the Nf - 2 
PQM model (left part) and those from the FRG (right part) for 
the chiral condensate (red and blue lines) and the Polyakov-loop 
expectation value (full lines) iflOl . The curves marked in blue dis- 
play pure QM results. 



field result for ^2lfTT1l. i.e. 

Q = Qmf + • 



(16) 



Diagrammatically this amounts to a summation of all ring 
diagrams (Fig. [3]) and represent mesonic contributions to 



^^^^^^^^ 



+ 




Fig. 3. Next- to-leading-order 1 /A^^ -corrections to the thermody- 
namic potential in the (P)NJL model. The wiggly lines denote the 
four-quark contact interaction of the model. These ring diagrams 
can be summed to all orders. 



Q in leading order: 



Qm 

M 



(17) 



In imaginary time Qm is represented by a Matsubara sum 
involving the mesonic polarization function TIm 



Q 



M 



r- 



2G77M(/^„q)) (18) 



which can be Wick-rotated and leads to 



'M 



■J (2^0^! V 



(1 + 2nB{(jj))^M 



(19) 



where is the finite temperature bosonic occupation prob- 
ability and (pM can be interpreted as an 'in-medium' phase 
shift in the Beth-Uhlenbeck sense. Explicitly one has 



^ 1 ^^ l-2G77M(^-/y/,q) 
2i ^ 1 - 2GnM(oj + IT], q) 



(20) 



The resulting mesonic contributions to the pressure in the 
three-ffavor PNJL model are displayed in Fig. [4] |[T2l . As 
mentioned above, the thermal quark- and diquark excita- 
tions are supressed by the Polyakov loop in the confined 
phase and the pressure comes entirely from mesons, as it 
should be. Up to temperatures of about 150 MeV it agrees 
well with the model-independent results from chiral per- 
turbation theory 1 13 |. When quarks and gluons take over 
above 200 MeV, the mesonic contributions become supressed, 
since the meson masses grow rapidly, reaching the ther- 
mal value of 27iT around 300 MeV. Hence they become 
strongly Boltzmann- supressed and the 'Hagedom singu- 
larity' is avoided. 

Now is is easy to see why the chiral quark-hadron tran- 
sition is softened when including ffuctuations beyond the 
mean ffeld. The chiral condensate is given by the derivative 
of the pressure (or Q) with respect to the bare quark mass: 
(qq) = dQ/dniq. Hence in the confined phase 



This softening can be understood physically by consider- 
ing the next- to-leading-order 1 /A^^ -corrections to the mean 
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Fig. 4. Various contributions to the normalized pressure PjT^ in 
the three-flavor PNJL model |12|. The dashed-dotted line indi- 
cates the mean-field result, while the dashed line displays the 
purely mesonic contributions (mostly pions). Up to temperatures 
of ~ 150 MeV these agree well with the results from chiral per- 
turbation theory L13J (dotted line). 



which is non- vanishing. Since the Polyakov loop dynami- 
cally couples to the quaks its softening can be understood 
by the same token. 



3 'Quarkyonic' Phase at low T and large yu? 

Based on large- A^c QCD it has been argued that there should 
exist a phase of strong-interaction matter in which the chi- 
ral and the deconfined transition split apart, since the de- 
confinement transition temperature becomes independent 
of the quark chemical potential 1 14 |. Within a narrow win- 
dow in II with a width of the order of ~ 1 IN^ there should 
then be a rapid transition to a dense phase of confined 
hadrons (mostly baryons at low T) in which chiral sym- 
metry is restored leading to parity doubled color singlet 
excitations. Schematically this leads to the following phase 
diagram: 



Within a PNJL model for Nc ^ oo and a fixed number 
of flavors this picture can be reproduced in the mean-field 
aproximation 1 15|. At A^^ = 3 the situation changes, how- 
ever. Fig. [6] showns a PQM calculation of the phase dia- 
gram for three colors and three flavors |6|. In the physical 
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Fig. 6. Mean-field results for the phase diagram in the Nf - 3 
PQM model |6|. The deconfinement boundary is determined by 
the inflection point in the Polyakov-loop expectation value i. 



word the (pseudo)critical transition lines curve downwards 
with increasing ji and the triple point evolves into the CEP. 
As one can see, there is a large region in which there still is 
confinement (in the statistical sence) but chiral symmetry 
is restored which thus qualifies as the quarkyonic phase. 
The size of this phase, however, crucially depends on the 
Polyakov-loop dynamics. 

Let us consider the polynomial ansatz for the Polyakov- 
loop potential for example 



(22) 



with 



b2{T) = flo + «i(ro/r) + a2(To/Tf + a^in/Tf . (23) 



T 



De-Confined 



Mesonic 



Quarkyonic 



Mb 



Fig. 5. Sketch of the large-A^c QCD phase diagram as it emerges 
from the arguments of Ref. 1 14 j. All phase boundaries are of first 
order with a 'triple point' where all first order lines meet. As a 
function of j^b the transition happens at the vacuum baryon mass 
Mb ~ NcAqcd- 



Originally the parmeter To was adjusted to the first-order 
deconfinement transition temperature ofT = 270 MeV in 
'pure gauge' theory, i.e. without dynamical quarks. This, 
however, does not take into account that with dynamical 
quarks Tq acquires a Nf - and fi dependence. Based on the 
one-loop running of the QCD yS-function this dependence 
has been estimated in Ref. Iil6j| as 



To(Nf,/i) = TrCxp{-l/aob(Nf,/j)) (24) 



with 



KA^/,yu)= ^(11A^,-2A^^)■ 



16A^ 



(25) 



where = a(A) is the running gauge coupling at some 
t/V-scale and Tj = 1.777 GeV is fixed at the scale of the 
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T-meson to reproduce Tq = 270 MeV at A^/ = with the 
corresponding value ao = 0304. 

At = this leads to the Polyakov-loop Tq given in 
Tab. 1 . Of importance for the size of the quarkyonic phase 
is the yu-dependence of Tq. When using the estimate in Eq. 



Nf 1 
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To [MeV] 
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Table 1. T and ytz dependence of Tq 



( [24| ) one obtains in a PQM calculation instead of Fig.[6]the 
result displayed in the upper part of Fig. [7] The region of 
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Fig. 7. Upper part: mean-field results for the phase diagram in the 
Nf = 3 PQM model |6| with To(Nf,/d). Lower part: FRG result 
for the Nf = 2 PQM model L17J . 



the quarkyonic phase has considerably shrunk. This situa- 
tion is even more dramatic if, instead of a mean-field calcu- 
lation, one includes quantum fluctuations via the FRG 1 17 1 
(lower part of Fig. [7]). In this case the chiral and the de- 
confinement transition essentially coincide and hence the 
region of quarkyonic matter has collapsed to nearly zero. 

There is a more heuristic argument, why the bound- 
aries for the chiral- and deconfinement transitions cannot 
be to far appart ifTSl . This is based on the statistical model 
of a hadron resonance gas which is extremely successful in 



describing the chemical freeze out of hadrons in relativis- 
tic heavy-ion collisions over a very wide range of beam 
energiesi 1 1. The thermodynamics of the resonance gas is 
that of a free gas and therefore the EoS can be straight- 
forwardly evaluated. Results for the entropy density s and 
the baryon number density n, normalized to the free quark- 
gluon gas values are displayed in the upper part of Fig. [8] 




Baryon Chemical Potential [MeV] 

Fig. 8. Upper part: The normalized entropy- and baryon number 
density of the hadron resonance gas as a function of T and 1 18J. 
The red and blue bands indicate the regions, where both quanti- 
ties increase from 0.3-0.8. The chemical freeze-out data are taken 
from II 191201 . Lower part: boundaries of the regions in which the 
chiral condensate and the Polyakov loop change between 0.4 and 
0.6. 



The idea is now to take the jd dependence of the Tq param- 
eter in the Polyakov-loop potential tl directly from the ji 
dependence of the empirical freeze-out line. With 

To(ji)/To(0) = 1 - (bTo)(jiB/Tof (26) 

and fixing Tq at 200 MeV this yields bTo = 2.78 x lO'^ 
MeV, in quite good agreement with the estimates of Ref. 
1 16 ] which give To = 187 MeV and bTo = 2.1 x 10"^ MeV. 
The mean-field PNJL calculation with this parametrization 
of Toijd) reproduces the entropy- and number density val- 
ues of the resoncance gas, shown in the upper part of Fig. 
[8] very well. The T and yu dependence of the chiral con- 
densate and the Polyakov loop can then be inferred from 
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the calculation. The results are shown in the lower part of 
Fig.[8]and again indicate that the chiral- and deconfinement 
transition line are close in the whole (r,yu)-plane, leaving 
little room for a quarkyonic phase. 



4 Inhomogeneous Phases of QCD Matter 

In the conventional picture which emerges from model cal- 
culations of the PNJL or PQM type, the chiral transition in 
the (r,yu)-plane is a smooth cross over at low yu and large 
T (see Sect. [2]). It becomes a first-order phase transition at 
large yu and small T from the chirally broken to the chirally 
restored phase. Usually it is assumed that both phases are 
homogeneous. From studies of dimensional fermionic 
theories, such as the Gross-Neuveu, the NIL or the t'Hooft 
model it is known analytically that in some regions of the 
(r,yu)-plane, inhomogeneous phases are preferred in the 
N ^ 00 limit | 21 1. When applying these results in a Ginz- 
burg-Landau analysis near the CEP (which is a second- 
order phase transition) it is found that the homogeneous 
phase in three space dimensions is unstable against one- 
dimensional spatial oscillations of the chiral order param- 
eter |22J. 



4.1 Phase diagram in 1+1 dimensional theories 

Let us consider the Gross-Neveu (GN) model with N de- 
grees of freedom as an example. With strong-interaction 
physics it shares some essential features, such as asmp- 
totic freedom, dimensional transmutation and spontaneous 
chiral symmetry breaking in the vacuum. The GN model is 
specified by the following Lagrangian: 



£gn = Yj ^^'^ - ^0) ^^'^ + y Z 



(27) 



where i// is an A/^-component fermion field, the coupling 
constant of the four-fermion interaction and mo the bare 
fermion mass. The model posesses a global U(N) symme- 
try and a Z2 chiral symmerty: 



(28) 



As usual, mean field theory becomes exact in the N 

00, Ng^ = const limit and the resulting self consistent Hartree 

equations for given T and ji read: 



{-iysd, + /M(z)) ij/a = eailfa 
M-mo = -Ng^ ^ n^ii/a^a • 



(29) 



Here na(T,fi) denote the Fermi-Dirac occupation probabil- 
ities. The mean-field equations can be solved analytically 
with the z-dependent (constituent) mass function given by 
an elliptic Jakobi-sn function in the chiral limit, mo = 0: 



The elliptic modulus v varies continuously between zero 
and one and ^ is a scale, related to the maximum of M(z). 
For y = I one has M(z) = ^tanh(^x), i.e. a single soli- 
ton and for y ^ the shape becomes more and more si- 
nusoidal, albeit the amplitude also goes to zero. Thus the 
sn-function interpolates smoothly between soliton-like and 
sinusoidal shapes. 

The actual values for given T and ji are determined 
from the minimization of the grand potential 

Q(T,j^) = -TTrlog{s-')) + £dz M{zf (31) 

where S is the fermion propagator in the Hartree approx- 
imation and A the period of the spatial modulation. The 
resulting phase diagram in the chiral limit is shown in in 
the upper part of Fig. |9] Raising T for jd < 21 n one en- 




M{z) = ^^vq sn(^z|v) 



(30) 



Fig. 9. Upper part: phase diagram of the exactly solvable 1 + 1 
dimensional Gross-Neveu model in the large-A/^ limit | 21 1. A spa- 
tially ordered phase occurs for small temperatures and 11 > Ijn. 
The phase boundaries deliniate second-order transitions which 
meet at a tricritical point, the 'Lifschitz' point. Lower part: evolu- 
tion of the spatial modulations of the mass function M(z) as one 
moves along a line of constant (small) T . While at smaller \i the 
modulation is solitonic it becomes increasingly sinusoidal with 
increasing \i. 



counters a line of second order transitions from the chi- 
rally broken, M 9^ 0, to the chirally restored phase, M = 0, 
both of them spatially homogeneous. Cutting through the 
phase diagram at small constant T, one enters a region of 
inhomogeneous phases through a second-order transition, 
where the spatial modulation go from soliton-like to sinu- 
soidal as \i increases (lower part of Fig. [9]). At large enough 
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fi and moderate T one again enters the homogeneous chi- 
rally restored phase through a second order transition. All 
lines of second order transitions meet in a tricritical point, 
which is commonly referred to as the 'Lifschitz' point. 



As one goes away from the chiral limit (lower part of Fig. 



TQ| ) the Lifschitz point moves to higher ji and lower T and 
the inhomogeneous region shrinks^ The phase diagram in 



4.2 Inhomogeneous phases in the NJL model 

As mentioned above, QCD-inspired models of the NJL or 
QM type are currently used to asses details of the phase 
diagram at large chemical potentials, where ab-initio lat- 
tice methods fail. It is therefore interesting to ask whether 
spatially inhomogeneous regions also occur in three spatial 
dimensions and what their properties are. This question has 
been addressed recently for ID modulations of the chiral 
order parameter in 3D space (plates) |23|. 

The starting point of the analysis is the two-flavor NJL 
model with scalar coupling: 

Znjl = q {W -m,)q^ Gs ({qqf + {qiy'r'q)^] • (32) 

Considering phases with a spatially varying expectation 
values {qq(x)} = S(x) and (iqiy^T^q(x)j = Pa(x) and re- 
stricitng to the case where the direction of the vector Pa(x) 
is constant in flavor space such that Pi(x) = P2(x) = 
and Psix) = P(x), the following mean-field Lagragian is 
obtained: 

4,j^(x) = q [W -m, + 2Gs {S (x) + iy\^P{x))) q 

-Gs{S(xf + P(xf) . (33) 
In terms of the (complex) mass function 



Mix) = mq- 2Gs (S(x) + iP(x)) 



(34) 



the mean-field grand potential for ID modulations can be 
evaluated straightforwardly and one obtains 



i2(r,,).-^yr^in 



2 cosh 



IT 



|M(x) - mj- 
4G.V 



(35) 



where V is the volume of the Wigner-Seitz cell of the peri- 
odic condensate. Here it is assumed that the inhomogene- 
ity is in the z-direction and the perpendicular (x, z/)-plane is 
translationally invariant. The discrete parameter a involves 
the Eigenvalues of the mean-field Hamiltonian in the z- 
direction which can be obtained analytically for the case 
in which M(x) is reaQ It can be argued that a real order 
parameter is thermodynamically preferred, at least in the 
vicinity of a second-order transition and in the chiral limit 
1221. 

As can be seen from Fig. [TO] the homogeneous first- 
order chiral transition is completely covered by an inhomo- 
geneous region, bounded by second-order transition lines. 



' A complex mass function M(x) leads to so-called 'chiral spi- 
rals' for sinusoidal modulations |25| 
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Fig. 10. Upper part: phase diagram of the Nf - 2 NJL model 
in the chiral limit | 23 1. The shaded area denotes the region of ID 
plates in the z-direction. The (tricritical) CEP of the homogeneous 
phases coincides with the Lifschitz point. The dashed lines denote 
the spinodal lines of the homogeneous phases. Lower part: same 
as the upper part, but now including finite bare quark masses of 
mq^5 MeV and 10 MeV 1231 . 



the 3-1-1 dimensional NJL model looks very similar to the 
case of the GN model, except that the inhomogeneous 
region does not extent to large jd. This may be due to cut- 
off' eff'ects |24|. While the l-i-l dimensional GN model is 
renormalizable, this is not the case for the NJL model in 
three dimensions. Here the loop intergrals have to regular- 
ized and the ultraviolet cut-off' enters as an explicit param- 
eter into the calculation. 



4.3 Including Vector Interactions 

It is well known that extended (P)NJL models which in- 
clude vector interactions lead to substantial modifications 
of the location of the CEP in the phase diargram. Eor suflfi- 
ciently large vector-interaction strength, the CEP even dis- 
appears and the chiral transition becomes a cross over in 
the whole (r,yu) plane. In view of the above discussion is 



The same conclusions also hold in the QM model iBl 
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therefore important to ask, how spatial inhomogeneities of 
the chiral order parameter influence these findings. 

Including vector interactions the Nf = 2 NJL Lagrangian 
reads 



£njl = q {W - m,) q+Gs ({qqf + [qij' r"" q)^yGy{qy^ qf . 

(36) 

Varying Gy strongly aff'ects the location of the CEP when 
only homogeneous phases are considered (upper part of 



Fig. 11). 
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Fig. 11. Upper part: The homogeneous phase diagram including 
vector interactions of increasing strength (left to right). At suffi- 
ciently large coupling the CEP disappears. Lower part: same but 
allowing for ID inhomogeneities L29 J . 



One-dimensional spatial inhomogeneities can be treated 
in a similar way as for the NJL model with scalar inter- 
action only. An additional complication arises, however, 
from the fact that the vector interaction couples to the den- 
sity which is also spatially modulated, i.e. n = i^qy^q^ = 
n(x). Via the renormalized quark chemical potential 



this also induces a spatially varying chemical potential. 
This complicates the minimization of the grand potential 
considerably. In order to take advantage of the fact that, 
for ID modulations, the partition sum of the NJL model 
without vector interaction can be performed analytically 
for states in the z-direction, one has to replace density by 
its spatial average 



^(x) n = {n(x)) - const. 



(38) 



/i(x) = - 2Gvn(x) 



(37) 



which is only approximate, but can be justified rigourously 
in the vicinity of a second-order boundary to the restored 
phase and in particular near the Lifschitz point |29|. As a 
consequence // becomes constant as well and the problem 
reduces to the case discussed in the previous Sect., albeit 
with shifted chemical potential fi. Another interesting con- 
sequence of non-vanishing Gy is that the CEP and Lifs- 
chitz point no longer coincide. As shown in the lower part 
of Fig.[TT] increasing Gy leads to two eff'ects: (1) the CEP 
moves deeper into the inhomogeneous region and even- 
tually dissappears, (2) the inhomogeneous region grows . 
The latter is largely a trivial eff'ect of the shifted chemical 
potential. Plotting the phase diagram in terms of the aver- 
age density, instead of p., the transition lines do not depend 
on Gy at all. It is remarkable that, although the CEP is 
strongly eff'ected by the vector interaction, the location of 
the Lifschitz point is not (except for a more or less triv- 
ial yu-dependence). This speaks for the 'robustness' of the 
inhomoneneous region. 

The influence of gluonic degrees of freedom on the oc- 
curence of spatially modulated phases via the coupling to 
the Polyakov loop has also been considered in Ref. |29l. 
Preliminary results indicate that the general picture remains 
unchanged. 

It has long been known that ID long-range order is 
thermodynamically unstable II26I27L In the context of pion- 
condensed phases this has been reiterated by Baym et al. 
|28 1. The occurence of ID order is a mean-field or laigc-N 
artefact and any finite temperature destroys the long-range 
order. There remains however the possibilty of a quasi- 
ordered one-dimensional phase, with long-range correla- 
tions decaying algebraically in space [28]. In any case, 
two- and three-dimensional structures will be stable and 
it is worth exploring them ftSQl . 



5 Conlusions and Outlook 

Much progress has been made in recent years in the theo- 
retical exploration of the phase diagram of hadronic mat- 
ter. Largely based on calculations in QCD-inspired mod- 
els a rich structure has emerged, especially at large chemi- 
cal potentials and low temperatures. With the help of ERG 
methods is has become possible the include fluctuations 
on top of the mean-field results 1191101171 . These are able 
to properly describe the critical exponents at the CEP f9\, 
lead to a small critical region around the CEP and result in 
a softening of chiral and deconfinement transition as com- 
pared to mean-field predictions at small yu. The latter find- 
ing is consistent with state-of-the-art lattice simulations of 
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the quark-hadron transition at vanishing ju and has a sim- 
ple physical interpretation in terms of meson fluctuations 
below the pseudo-critical temperature |[T2l. 

An interesting question that has caught much attention 
recently is the existence of a quarkyonic phase which is 
confining but chirally restored. The arguments for such a 
phase are based on a large Nc analysis of the pressure at 
finite T and ji. For Nc = 3 there are substantial modifica- 
tions and it is not clear in what form the large-A^^ picture 
survives. Mean-field PNL and PQM calculations indicate 
that there is a region in the phase diagram in which quarky- 
onic matter could exist. The size of this region crucially 
depends of the parameters of the Polyakov-loop potential, 
in particular the yu-dependence of Tq. Present FRG calcu- 
lations indicate that the deconfinement and chiral transi- 
tion lines in the (T,jd) plane almost coincide 1 17 L leav- 
ing little room for a quarkyonic matter state. These find- 
ings are corroborated by a heuristic analysis based on the 
resonance-gas EoS and the experimental location of the 
chemical freeze-out points in the phase diagram flSl . 

An exciting new possibility is the occurence of spa- 
tial modulations of the chiral order parameter, leading to a 
ordered density profile. Exact dimensional results for 
QCD-like models in the large-A^c limit can be taken over to 
three space dimensions in the PNJL and the PQM model 
|23 1 and lead to plate-like structures. The inhomogeneous 
phases are bounded by second-order transition lines and 
lead to a Lifschitz point in the chirial limit. Without vector 
interactions, the Lifschitz point coincides with the CEP but 
both drift appart with increasing vector coupling 1291 . One 
can even achieve scenarios in which the CEP dissapears 
while the inhomogeneous phase region remains. It is very 
robust and its size is unaff'ected when considering the aver- 
age quark-number density rather that ji. One-dimensional 
phase are thermodynamically unstable although quasi one- 
dimensional structures ar still possible. Two- and three di- 
mensional spacial ordering is likely to occur and should be 
studied 1301 . It is well established theoretically that such 
phases are posible near the crust-liquid interface in a neu- 
tron star and it would be very exciting if they would also 
show up in the quark-hadron transition in its inner core. 
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